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Abstract 

The ground state phases of a one-dimensional SU(4) spin-orbital Hamiltonian in a generalized 
external field are studied on the basis of Bethe-ansatz solution. Introducing three Lande g factors 
for spin, orbital and their products in the SU(4) Zeeman term, we discuss systematically the various 
symmetry breaking. The magnetization versus external field are obtained by solving Bethe-ansatz 
equations numerically. The phase diagrams corresponding to distinct residual symmetries are given 
by means of both numerical and analytical methods. 
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I. INTRODUCTION 



There have been much interest in the study of spin models with orbital deEfeneracy 
due to experimental progress related to many transition-metal and 
rare-earth compounds such as LaMnOa , CeBg, and perovskite lattice, as in KCuFa 
Those systems involve orbital degree of freedom in addition to spin ones. Almost three 

10| had pointed out the possibility of orbital excitations in 



decades ago, Kugel and Khomskii 



these systems. As a model system, it exhibits some fascinating physical features which is lack 
without orbital degree of freedom. The isotropic case of spin system with orbital degeneracy 
was shown to have an enlarged SU(4) symmetry and one dimensional model is known to 
exactly solvable 0,0]. Materials related to spin-orbital systems in one dimension include 
quasi-one-dimensional tetrahis-dimethylamino-ethylene (TDAE)-C6o[12], artificial quantum 
dot arrays!^, and degenerate chains in Na2Ti2Sb20 and Na2V205 compounds |l^ lisj ]. 
It is therefore worthwhile to systematically study the features of one dimensional model. 
Theoretical studies P] has found the strong interplay of orbital and spin degrees of freedom 
in the excitations spectra. It has been noticed that the presence of orbital may results 
in various interesting magnetic properties. Applying a conventional magnetic field, the 
spin orbital chain with SU(4) symmetry is shown to reduce to a model with orbital SU(2) 
symmetry|S] in the ground state. Recently, we showed that the magnetization process 
becomes more complicated if taking account of the contribution from orbital sector We 
have explained that the competition between spin and orbital degree of freedom leads to an 
orbital anti-polarization phase. However, the external field we introduced in ref is not the 
most general one for SU(4) systems. 

In this paper, we study a SU(4) spin-orbital chain in the presence of a generalized external 
field on the basis of its Bethe ansatz solution. Our paper is organized as follows. In Sec. 
im we introduce the Bethe-ansatz solution and the Zeeman like term which is going to be 
added to the original SU(4) Hamiltonian. In Sec. IIIH we give some useful remarks on the 
quantum number configurations for the ground state in the presence of external field that is 
characterized by three parameters. We also simply demonstrate the thermodynamic limit of 
the Bethe-ansatz equation and briefly present the dress energy description of ground state 
in the presence of external field. In Sec. IIVI we study the magnetization properties of a 
Hamiltonian in the regime with one-parameter symmetry breaking. In Sec. we study 
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both magnetization and the phase diagram in the regimes with two-parameter symmetry 
breaking. Various phases and the quantum phase transitions are obtained by both numerical 
calculation and analytical formulation. Concerning to various phases we present detailed 
explanation in terms of group theory. Sec. I VII gives a brief summery. 



II. THE MODEL AND ITS SOLUTION 



We start from the following Hamiltonian 

where Sj and Tj denote respectively spin and orbital operators at site j, both are generators 
of SU(2) group characterizing the spin and orbital degree of freedom of outer shell electrons 
in some transitional metal oxides at the insulating regime. The coupling constant is set to 
unit for simplicity. It has been pointed that the above Hamiltonian possesses an enlarged 
SU(4) symmetry|l| rather than SU(2) xSU(2) symmetry. 

The four states that carry out the fundamental representation of SU(4) group is denoted 

by 

IT) = |l/2,l/2), IT) = 11/2,-1/2), 

li) = 1-1/2,1/2), II) = 1-1/2,-1/2) (2) 

These bases are labelled by the eigenvalues of and T^, i.e., \S^,T^). As the sm(4) 
Lie algebra is of rank 3, there exists third generator 2S^T^ which possesses simultaneous 
eigenvalue together with and T^. For convenience, we denote this new generator by 
hereafter. In the terminology of group theory, however, the quadruplet can also be 
labelled by the weight vectors which is defined by eigenvalues of Of, O2, 0| that constitute 
the Cartan subalgebra of su{4) Lie algebra. Here we adopt the Chevalley basis because the 
physical quantities can be conveniently expressed in this basis. 

The eigenvalues of S^, T^, as well as that of Of, Of, 0| are given in Table I, the relation 
between these two basis readsQl, 

S' = 01 + 20^2 + 01 
T' = 01 + 01, 

= 01-01 (3) 



TABLE I: The eigenvalue of S^, T^, and z-component of Of^,0|,0| for the four basis states 
(Eq. H. 

State Of Of Of 

It) 1/2 1/2 1/2 1/2 

IT) 1/2 -1/2 -1/2 -1/2 1/2 

li) -1/2 1/2 -1/2 -1/2 1/2 

II) -1/2 -1/2 1/2 -1/2 



The present model has been solved by Bethe-ansatz method ^jIUIj its energy spectrum 
is given by 

M 

Eo{M,M',M") = -J2yj^- (4) 

a=l ' 

where the A's are solutions of the following coupled transcendental equations 



M 



a'=l 

M' 

6=1 

M M' 

27rJb = y^^9-i/2il^b - Ag) + y^6'i(/ife - Hb') 

a=l b'=l 
M" 

c=l 

M' M" 

2itKc = ^9-1/2(1^0- IJ'b) +'^di{i^c- l^c'), (5) 
6=1 c'=l 

where 6a{x) = — 2 tan^^(x/a). The A,/i and u are rapidities related to the three generators 
of the Cartan subalgebra of the su{A) Lie algebra. The quantum numbers Jb, K^} specify 
a state in which there are N — M number of sites in the state ||), M — M' in ||), M' — M" 
in II), and M" in |J,). Hence the z-component of total spin, orbital and f/^ are obtained as 
Sl^ = N/2 - M\ T^^ = N/2 - M + M' - M", and U^l^ = N/2-M + M" . 

In present SU(4) model, a three-parameter external field (/ii,/i2,^3) can be introduced 
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to write out a most general Zeeman-like energy. 

3 

'Hzee = ^ hmO^. (6) 
m 

For more clear physics implication, we re-choose the parameters to write the effective mag- 
netization Aiz, 

M^=gsS!,, + g,T,l, + g^U,l,. (7) 

where gg, gt, gu are generalized Lande g factors for 5^, and t/^ respectively. Eq. ((7j) can 
be expressed in terms of the number of rapidities, 

N 

= —{gs + gt + gu)-M{gt + gu) 

-M\gs-gt)-M"{g,-g^), (8) 

Because the Zeeman-like term commutes with the SU(4) Hamiltonian (^, the energy 
spectrum in the presence of external field are simply related to the energy spectrum in the 
absence of external field, 

E{h, M, M', M") = Eo(M, M', M") - hM', (9) 

where Eo{M, M', M") is determined by eqs. (jSj). Obviously, the apphcation of the external 
field with different magnitude just brings about various level crossings. 
In terms of Oi, O2, O3, the magnetization ((Tj) becomes 

= {gs + gt + gu)0{ + 2gsO'^ + {g^ + gt- gu)Ol. 

(10) 

which breaks SU(4) symmetry down to various lower symmetries depending on the distinct 
regions in the parameter space. 

III. THE GROUND STATE CONFIGURATION 

Based on the Bethe ansatz solution of the model, we first give the quantum number 
description of the ground state, which is useful for numerical method. We also give the dress 
energy description for the ground state and propose the conditions to determine quantum 
phase transitions, which is useful for analytic study. 
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It has been known that the ground state of the Hamiltonian (^J) is a SU(4) singlet for 
the case of = An. The configuration of the quantum number for ground state {la, Jb, Kc} 
(a = 1,2, ...,3n; b = l,2,...,2n; c = l,2,...,n) are consecutive integers (or half integers) 
arranging symmetrically around the zero. In the presence of magnetic field, however, the 
Zeeman term brings about level crossings and the state with M = 3n, M' = 2n, M" = n is 
no longer the ground state. Therefore the numbers M, M', M" for the lowest energy state 
are related to the magnitude of the applied external field. 

In order to solve the Bethe ansatz equation numerically, we need to determine the possible 
configuration of quantum numbers for give values M, M' and M" . The property of Young 
tableau requires that Max(M) = 3A^/4, Max(M') = N/2, Max(M") = and - M > 
M — M' > M' — M" > M" for a given A^. Then one is able to analyze the change of energy 
level for each state, which determines the true ground state for a given external field. One 
can also calculates the magnetization by Eq. (jHJ. 

In the thermodynamic limit, the energy (jH]) is expressed in terms of densities of the 
rapidities. 



EIN =--{gs + gt + gu) 




J-Xo 




-Mo 




(11) 
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These densities satisfy the following coupled integral equations 



a 



(A) = iri/2(A) - K^{\ - \')a{X')dX' 

/MO 
Ki/2{X - fi)uj{fi)dfi, 
Mo 



Ki/2{fi - \)a{\)d\ 

-An 

Mo 

/I'O 
Ki/2{n - u)r{u)du, 
-UQ 

/MO 
Ki/2{u - fi)uj{fi)dfi 
-Md 



MO 



Mo 

I'O 

iri(z/- z/')r(z/')c/z/'. (12) 

-1^0 



where Kn{x) = n ^n/{n'^ + x^) and Ao,/Uo and z/q are determined by 



M' 





a(A) 


>/-Ao 




/•Mo 




''-MO 




/•I'O 









(13) 



AT 

It is more convenient to introduce dress energy 16]. The iteration of eg . (1121) gives rise to 
e{\) = -2TTKi/2iX) + igt + 9u)h 

"A,, 



- / " /^i(A- A')£(A')rfA' 

/Mo 
Ki/2{\ - ^i)C{fi)dfi. 
-Ul) 



-Mo 

/Ao 
Ki/2{iJ, - X)e{X)dX 
Ao 

Mo 

Ki{fi - fi')C{fi')dfi' 

MO 

Ky2ifi - iy)^{iy)diy 

/•MO 

^{i^) = {gt - gu)h + / Kii2{v - fi)({fi)dfi 



MO 



/I'O 
K^iu-u')aly')du'. (14) 



where e, (, and ^ are the dress energies in A, /i and u sectors respectively. It is worthwhile 
to point out that the dress energy is also the thermal potentials at zero temperature, i.e., 
exp(£:/T) = p^/p, exp((^/T) = a^/cr, and exp(^/T) = uj^/uo. The thermal Bethe-ansatz 
equation in the zero-temperature limit turns to eq. (|14j) . In terms of dress energy, the energy 
(fTT|) is simplified to 

E/N = + 9t + gu)+ r K„^{\)e{\)d\. (15) 

Apparently, the ground state is a quasi-Dirac sea where the states of negative dress energy, 
e{X) < 0, Q{p) < 0, ^{v) = 0, are fully occupied. The Fermi points of the three rapidities 
are determined by 

£(Ao) = 0, C(/^o) = 0, e(z/o) = 0. (16) 

The system will be magnetized if the applied field enhance the dress energy, because it makes 
the corresponding Fermi points decline. The quantum phase transition occurs when any of 
the Fermi points shrinks to zero. As a result, the critical values of the external field are 
solved by 

£(0)1,=,, = 0, cmh=K = 0, m\h=h'^ = (17) 

These conditions together with eq. ()14|l enable one to calculate those critical values. 

IV. REGIMES WITH ONE-PARAMETER SYMMETRY BREAKING 

The application of external field makes the SU(4) symmetry break down to various 
regimes with different residual symmetry. In this section, we shall discuss the simplest 
cases of single parameter hierarchy. There are three special directions in the weight space of 
SU(4). If the external field is supplied along those directions, i.e., either hi, h2 and in eq. 
(jni) is not vanished, a partial breaking of a SU(2) to U(l) will take place. Let us consider 
them respectively. 
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A. Residual SU(3)xU(l) symmetry 



If Qs = 0,gt = —Qu > 0, the Zeeman interaction pO|) becomes 



= 2gtOl 

= gtM' - 2gtM". (18) 

The occurrence of operator Of makes the Hamiltonian noncommutable with Of. Thus a 
SU(2) subgroup generated by 0|, is broken down to U(l). Analyzing the level 

crossing from Eq. (jU)), we shown the magnetization curve in Fig. ^ 

Because the terms gt+Qu and gg—gt in the first two equations of eqs. (fT^ are non-positive, 
the external field can not enhance the two dress energies e{X) and C(/^)) the Fermi points in 
both A and /i sectors are fixed. On the contrary, the Zeeman term has a positive contribution 
to and its two Fermi points will decline when the external field increases. Although it 
is a SU(4) singlet labelled by Young tableau [n^] in the absence of external field, the ground 
state possess a residual SU(3)xU(l) symmetry in the presence of the aforementioned one 
parameter external field at small magnitude, which corresponds to phase IV labelled by 
four-row Young tableau. In this regime are there still three type of rapidities that solve the 
Bethe-ansatz equation. The U(l) is generated by Of, while the SU(3) is generated by the 
following eight operators 

1 1 



01 = 




01 


= -(S' 


ot = 


ot 


= S+T 


or = 






= S'T 




3^1+ n, 




-2 = S' 



(19) 



There exists a critical field when those two Fermi points shrink to zero, the rapidity 
disappears in the Bethe-ansatz equation. Thus a quantum phase transition occurs at the 
critical field which separate two phases, we call phase IV and phase III. 

The magnetization process can be clearly illustrated by the evolution of Young tableau. 



which shows the evolution of Young tableau in magnetization process, i.e., from SU(4) 
singlet to SU(3)xU(l) states and then to a SU(3) singlet. Physically, we have M.^ jN = at 
zero external field due to there are a quart of the total sites being respectively in the states 
It), It), It) and It). Turning on the external field leads to the spin-orbital flipping, |t) |t), 
It) It) and It) It) I which result in non- vanishing magnetization. The SU(3)xU(l) 
symmetry makes the above three flipping processes favor to occur simultaneously. When 
the external field succeeds a critical value, it goes into phase III where all the states |t) have 
been flipped over. In this phase the 2;-component of total spin and total orbital keep positive 
constant S^/N = T^/N = 1/6 while that of keeps a negative constant U^/N = —1/6. 
Consequently, the magnetization reaches a saturate value A4^/N = 2/3 and the ground 
state becomes the SU(3) singlet regardless of the magnitude of the external field in phase 
III. 

B. Residual SU(2)xSU(2) symmetry 

Applying the external field along the direction of the second simple root of su{4) Lie 
algebra, we will have a symmetry breaking from SU(4) to SU(2)xU(l)xSU(2) for ground 
state. This is realized by the choice of Lande g factors = 0,gs = —Qt, which makes the 
magnetization to be 

= gM - 2gM' + 9sM". (20) 
Those two SU(2) are generated respectively by 

|i(T^ + f/^), (1 + S^)T±|, 

|i(T^_t/^), (i_^^)T±|. (21) 

As gs — gt in the second equation of (fT^ is positive but both gt + g^ and gt — gu in the 
first and third equation are negative, the external field makes the Fermi points in sector 
to shrink. The critical value of the external field when quantum phase transition occurs is 
determined by C(0)Uc = 0- This critical point separates two different phases, we call phase 
IV and phase II. 
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The magnetization curve is shown in Fig. It can also be illustrated by the evolution 
of the Young tableau, 



;;;b-b;;;;;;b 



The spin-orbital flipping process caused by the applied external field has several character- 
istics. During the flipping process, ||) and ||) flips simultaneously into |t) and |t) pairs, 
which makes four-row Young tableau reduce to the two-row Young tableau when across 
the critical field. Apparently, the eigenvalues of both and do not change during the 
magnetization process. Only contributes to the magnetization Ai^. The total spin is 
completely polarized (i.e., Ai^ is saturated) once the phase IV transits to phase II. The 
phase that Yamashita et a/ jsj discussed agrees with this special case. 



C. Residual U(l)xSU(3) symmetry 

If (7s = and gt = Qu, the magnetization becomes 

= 2gtOl 

= gt{N -2M + M'). (22) 

which implies that the external field was applied along the first simple root of s-u(4) Lie 
algebra. This gives rise to a symmetry breaking down to U(l)xSU(3) for ground state. For 
the sake of saving space, we omitted the operators that generate those symmetry. Because 
this parameter choice implies that gt + in A sector positive but the g factor terms in both 
and V sectors are non-positive, the quantum phase transition is only related to A sector. 
The critical value is determined by £:(0)|/i^ = 0. This critical point separates the system into 
two phases, phase IV and phase I. The magnetization process is shown in Fig^ 
It is helpful to illustrated this process by the evolution of Young tableau, 



i 




























i 





























From Fig^ we see that the spin, orbital as well as are polarized simultaneously versus 
external field. The above Young tableau indicates that the states |t), ||), ||) change to the 
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state II) simultaneously because they carries out a SU(3) representation and the system 
possess SU(3) symmetry. After the system is fully polarized, the magnetization reach the 
maximum value M.^ jN = 1/2. Then the residual symmetry of the ground state is only of 
U(l). 

V. REGIMES WITH TWO-PARAMETER SYMMETRY BREAKING 

In the previous section, we discussed the simplest case where merely one SU(2) subgroup 
symmetry is broken. In the following, we will consider the regimes with two-parameter 
symmetry breaking, which involves more SU(2) subgroup. 

A. Residual U(l)x U(l) x SU(2) symmetry 

Under the restriction Qu = Qs + Qt, the magnetization (jH)) becomes, 

M' = 2{g, + gt)Ol + 2g,Ol (23) 

which indicates that the residual symmetry of ground state is U(l)x U(l) x SU(2), in which 
the SU(2) is generated by Of = (T" - U')/2 and Of = T±(l/2 - S'). 

The magnetization curves for different gt is plotted in Fig. |21 and the phase diagram 
in terms of gt/gs versus h is given in Fig. IHl On the one hand, because the eigenvalue 
of equals that of for both states |t) and |t), the flipping process occurred in phase 
II contributes to both magnetization of and equivalently. On the other hand, the 
flipping from |J,) and ||) occurs simultaneously in the phase IV due to the SU(2) symmetry. 
As a results, the magnetization of and are expected to be the same in the whole 
process, which can be seen from our numerical calculation in Fig. |2l 

For gt/gs < 0.5, there exists three distinct phases, denoted by IV, II and I respectively 
according to the number of the rows of Young tableau. The SU(2) symmetry makes the 
states II) and ||) flip simultaneously when the external field increases. This makes the four- 
row Young tableau turn to two-row Young tableau directly, hence the phase III labelled by 
three-row Young tableau will not take place. The boundary between phase IV and phase II 
is determined from C(0) = and ^(0) = together, 

9t/9s = l + ^K^/2{0)*eiO), (24) 
12 



where e can be computed from eqs. ()14|) numerically. For sufficient large external field, all 
5, T and U are frozen to the z-direction, which brings about the occurrence of phase I. The 
boundary between phase I and phase II is determined by e{Q) = 0, i.e., 

= i - i (25) 

This can also be derived from the competition between the states related to the Young 
tableau [A^ - 1, 1] and [N]. 

The asymptotic behavior at large h is gt/gs = —1/2; which implies that the phase I will 
never occur as long as gt/gs < —1/2- The magnetization process in the region —1/2 < 
gt/gs < 1/2 can be illustrated by the following evolution of Young tableau, 

i 
i 

The boundary between phase IV and phase I is determined by e{0) = 0, C(0) = and ,^(0) 
together, 

9t/9s = 1^-1- (26) 

The common solution of Eqs. ()24II26|) gives h = 1 and gt/gs = 1/2, which is a three- phase- 
coexist point. 

B. Residual U(l)x SU(2) x U(l) symmetry 

If the external field along the direction of the second simple root is quenched but those 
along the other directions are kept, we will have symmetry breaking down to U(l)x SU(2) 
X U(l). Such kind of symmetry breaking is caused by Zeeman term of the following mag- 
netization 

M' = {gt + gu)Oi + {gt-gu)Os. (27) 

The magnetization curves with different gu is plotted in Fig. 01 and the phase diagram 
in terms of gu/gt versus h is given in Fig. El For gu/ gt < 1/2, there exists three phases 
denoted by IV, III and I respectively. 



i:::::i-i:::::i:P-'*a:::P-'*°-° 
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The boundary between phases IV and III is determined from ^ = 0, which can be solved 
numerically. For sufficient large external field, the magnetization is saturated reaching phase 
I. This phase transition occurs at 

9u/gt = \-\ (28) 

Thus if Qui Qt < the phase I will never occur regardless of the magnitude of external 

field. So in the region —1/2 < gu/ gt < 1/2, the magnetization process can be illustrated by 
Young tableau 



In phase III, the length of the second row and that of the third row in their corresponding 
Young tableau is always equal due to the SU(2) symmetry. Thus the probability of pure 
spin-fiipping |J,) — > |t) and pure orbital-fiipping ||) ||) is the same. Additionally, the 
process ||) |t) contributes the same for the magnetization of and T^. These properties 
result in the same magnetization curves of and shown in Fig. 0] Flipping over the 
state I J,) which has positive eigenvalue of brings about a negative magnetization of 
in phase IV. 

The phase I will never occur for gu/gt < ^1/2, as is similar to the case of SU(3) xU(l). 
Actually, it recovers the case of residual SU(3)xU(l) symmetry at gu/gt = —1, 

If gu/gt > 1/2, we can see from Fig. Elthat the phase IV transits into phase I directly 
along with the increase of external field. The boundary which separates these two phase is 
determined by 

9u/gt = ^ - 1- (29) 

Obviously, the point at gu/gt = 0.5 and /i = 2 is a three-phase-coexist point. The mag- 
netization properties in the region of gu/gt > 1/2 is similar to the case of U(l)xSU(3), 
particulary, a larger U(l)xSU(3) symmetry remains for gu/gt = 1- 
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C. Residual SU(2)x U(l) x U(l) symmetry 

The third case of two-parameter symmetry breaking is produced by the Zeeman term 
with the restriction to the the Lande g factor = —Qs — Qt- The magnetization now reads 

M' = 2g,Ol + 2{g, + gt)Ol. (30) 

which breaks the SU(4) symmetry down to SU(2)x U(l) x U(l). 

The magnetization curves for different gt/gs is plotted in Fig. IHl and phase diagram in 
terms of gt/gs versus h is given in Fig. [7| In the region of gt/gs > —1/2, there exist three 
phases denoted by IV, III, and II respectively. The final state is characterized by two rows 
Young tableau at sufficient large external field due to the SU(2) symmetry. The boundary 
separates phase III and II is determined by C(0) = 0, which reads 

9t/9s = 2 - — • (31) 

Obviously if gt/gs > 1/2, the phase II will never occur regardless the magnitude of external 
field, and the magnetization process is similar to the case of SU(3) xU(l). 

In the region —1/2 < gt/gs < 1/2, the magnetization process can be illustrated by the 
following Young tableau. 



Fig. ini shows that that the magnetization process of S^, and are quite different. In 
phase IV, the flippings from |J,) to other three states possess positive contribution to and 
T^, but negative contribution to U^. In phase III, undergoes polarization continually, 
while undergoes polarization but undergoes anti-polarization. Since in the final 
phase, both eigenvalues are zero and the spin magnetization are saturated, this phase will 
not change for any further increase of the external field. 

The point of gt/gs = 0.5, = ln2 is a three phase coexistence point. It can be seen that 
there exist two phases merely for gu/gt < —1/2, and phase IV transits into phase II directly 
. Actually, the residual symmetry of SU(2)xU(l)xSU(2) is restored when gt/gs = — 1- 
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VI. BRIEF SUMMARY 



In above, we studied the magnetization properties of a SU(4) spin-orbital chain in the 
presence of a generahzed external field that has three parameters due to the Cartan subalge- 
bra of su{4:) Lie algebra has three generators. Those three parameters are re-chosen so that 
to relate them with the spin S, orbital T and their product U. We called them three Lande g 
factors. Then all possible symmetry breaking and the corresponding magnetization process 
induced by that external field are studied respectively. The ground state phase diagram 
caused by the competition of quantum fluctuation and Zeeman-like effect are studied by 
solving the Bethe-ansatz equations numerically. The phase transition boundaries derived 
by studying the dress energy equations analytically. The features of various phases and 
transitions between them are explained in detailed in terms of group theory analysis. 

This work is supported by trans-century projects, Cheung Kong projects of China Edu- 
cation Ministry, and NSFC No. 10225419, No. 90103022. 
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FIG. 2: The magnetization M.^,Sz,Tz,Uz of the system with gs = 2,gu = gs + gt, and 
gt = 2.0, 1.5, 1.0, 0.5, 0.0, -0.5, -1.0. 
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FIG. 3: The phase diagram of ^t/gis versus h with gu = gs + 9t and residual symmetry U{1) x 
U{1) X SU{2) 
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FIG. 4: The magnetization Ai^,Sz,Tz,Uz of the system with gg = 0,gt = 2, and 
gu = 2.0, 1.5, 1.0, 0.5, 0.0, -0.5, -1.0. 
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FIG. 5: The phase diagram of gu/gt versus h with gs = and residual symmetry U{1) x SU{2) x 
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FIG. 6: The magnetization ^ Sz^T^^Uz of the system with Qs = 2,^^ = —Qs — gt-i and 
gt = 2.0, 1.5, 1.0, 0.5, 0.0, -0.5, -1.0. 



19 



0.5 



-1, 



- 




m 


1 

- 




\ 

\ 










n 




IV 




1 



external field: h 



FIG. 7: The phase diagram of gt/gs versus h with g^ = —gs — gt- and residual symmetry 
SU{2) X U{1) X U{1) 



20 



